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Abstract 

Background: The parasitic nematodes Ascaris lumbricoides and A suum are of great public health and economic 
significance, and the two taxa were proposed to represent a single species. miRNAs are known with functions of 
gene regulations at post-transcriptional level. 

Results: We herein compared the miRNA profiles of A lumbricoides and A suum female adults by Solexa deep 
sequencing combined with bioinformatics analysis and stem-loop real-time PCR. Using the A suum genome as the 
reference genome, we obtained 171 and 494 miRNA candidates from A lumbricoides and A suum, respectively. 
Among which, 74 miRNAs were shared between the two taxa, 97 and 420 miRNAs were A lumbricoides and A 
suum specific. Target and function prediction revealed a significant set of targets which are related to ovarian 
message protein, vitellogenin and chondroitin proteoglycan of the two nematodes. Enrichment analysis revealed 
that the percentages of most predicted functions of the miRNA targets were similar, with some taxon specific or 
taxon enhanced functions, such as different target numbers, specific functions (NADH dehydrogenase and electron 
carrier functions), etc. 

Conclusions: This study characterized comparatively the miRNAs of adult A lumbricoides and A suum, and the 
findings provide additional evidence that A lumbricoides and A suum represent a single species. Due to the fast 
evolution nature of miRNAs and the different parasitic living conditions of humans and pigs, the phenomenon 
above might indicate a fast evolution of miRNAs of Ascaris in humans and pigs. 
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Background 

Roundworms Ascaris lumbricoides and A suum are 
important parasites of human and pig health and socio- 
economic significance, with global distribution [1,2]. As- 
caris infects approximately 1.2 billion people globally 
and has been associated with intestinal pathology, re- 
spiratory symptoms and malnutrition in children from 
endemic areas [3]. It was reported that humans and pigs 
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can both be infected with the two nematodes [2]. There 
are very limited nucleotide differences in the sequences 
of internal transcribed spacers (ITS) of ribosomal DNA 
between the two ascarid species [4]. A. lumbricoides and 
A. suum are now proposed to represent a single species, 
and A, suum is considered a synonym of A. lumbricoides 
[5]. The genome and transcripts of A, suum were re- 
cently available, which provided valuable resources for 
better understanding and further studies of the biology 
of the Ascaris parasites [6-8]. 

MicroRNAs (miRNAs) are non-coding small RNA of 
18-24 nt in length. They are considered as key regula- 
tors for gene expression at the post-transcriptional level 
[9-12]. Due to their key regulating functions in growth, 
metabolism, development and cell differentiation and 
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their ability to respond to environmental and develop- 
mental signals, miRNA is essential for the complex life cy- 
cles of human and animal parasites [13-22]. miRNAs are 
also important for pathogen-host interactions [23-26] . Re- 
cently studies have indicated that miRNAs may represent 
potential biomarkers for various biomedical problems, 
such as the differentiation of different tumors and muscu- 
lar tissues [27-29], and they may provide biomarkers for 
the characterization of different genotypes of Toxoplasma 
gondii [30]. 

There might be differences in gene expression and regu- 
lation given the different living environments of A. lum- 
bricoides and A. suum, although they are considered to 
represent the same species. Therefore, it would be inter- 
esting to characterize the expression profiles of miRNAs 
in the two taxa. Herein the objective of the present study 
was to examine and compare the miRNA profiles of A. 
lumbricoides and A. suum using an integrative approach 
combining Solexa deep sequencing combined with bio- 
informatics analysis and stem-loop real-time PCR. 

Methods 

Ethics statement 

The present study was approved by the Ethics Committee 
of Lanzhou Veterinary Research Institute, Chinese Acad- 
emy of Agricultural Sciences (Approval No. LVRIAEC 
2011-006), and the A. lumbricoides and A. suum samples 
were collected strictly according to the requirements of 
the Ethics Procedures and Guidelines of the People's 
Republic of China. 

Nematodes 

Female adults of A. suum roundworms were obtained 
from slaughtered pigs in Shenzhen, China. Female adults 
of A. lumbricoides were obtained from a patient with as- 
cariasis after being treated with piper azine in Zhanjiang, 
Guangdong Province, China. Worms were incubated in 
physiological saline for 3 h at 37°C and then washed 3 
times to get rid of contamination from hosts. Female 
adults were identified morphologically and their identity 
was further ascertained by PCR amplification and se- 
quencing of the first internal transcribed spacer (ITS-1) 
of rDNA [4]. The nematodes were stored at -80°C 
until use. 

Total RNA isolation, small RNA preparation and deep 
sequencing 

Total RNA was prepared from a whole single adult of A. 
lumbricoides and A. suum respectively, using Trizol re- 
agent (Invitrogen) according to the manufacturers proto- 
col. Ten [ig total RNA and Novex 15% TBE-Urea gel 
(Invitrogen) were used for small RNA isolation. The RNA 
fragments of 18-30 bases long were added with 5' and 3' 
adaptors (Illumina) to the both ends, reverse transcripted 



with a RT-PCR kit (Invitrogen), and sequenced employing 
a Solexa sequencer at Huada Genomics Institute Co. Ltd, 
China. 

Computational analysis 

Adaptors, low quality sequences and reads smaller than 
18 nt were removed from the raw sequencing data. The 
reads were then searched against GenBank and Rfam 
database (http://www.sanger.ac.uk/software/Rfam) to iden- 
tify non-coding RNA. The remaining reads were mapped 
onto the A. suum genome by using the SOAP [31] with the 
sequences of pre-miRNA meeting the three criteria: 1) 
there was a standard stem-loop structure of pre-miRNA; 
2) mature miRNAs were present in one arm instead of 
the loop of hairpin precursors; and 3) the free energy 
hybridization was lower than -18 kcal/mol. The identified 
miRNA candidates were searched against the known miR- 
NAs of A. suum deposited in the Sanger miRBase with 
Blast [32]. Unmatched miRNA candidates were marked as 
novel miRNA. Targets of miRNA candidates were predi- 
cated with RNAhybrid [33] with the following extra pa- 
rameters: A) the AA G was set as lower than -25 kcal/mol; 
B) P-value was set as < 0.01. The Gene Ontology database 
(GO, http://www.geneontology.org/), Blast and Interpros- 
can (http://www.ebi.ac.uk/Tools/pfa/iprscan/) were used 
for prediction of functions of predicted targets. 

Analysis of novel miRNA expression 

The representative novel miRNA in the two nematodes 
were certified using stem-loop real-time reverse tran- 
scription polymerase chain reaction (RT-PCR) with 
SYBR Green [34]. The house keeping gene p-actin of A. 
suum (GenBank accession no. BI594141) was used as the 
endogenous control with primers as follows: forward pri- 
mer (5 -CTCGAAACAAGAATACGATG-3 ') and reverse 
primer (5 -ACATGTGCCGTTGTATGATG-3 '). Primers 
were synthesized by Shenggong Co, Ltd., China. The cycle 
conditions were as follows: 94°C for 30s, 52°C for 30s, 
72°C for 30s and finally with a single extension at 72°C for 
10 min. The quantification of each miRNA relative to (3- 
actin gene was calculated using the equation: N = 2" ACt , 
ACt = Ct miRN A-Ct acin [35] . 



Table 1 Common and taxon-specific reads of Ascaris 
lumbricoides and Ascaris suum 



Classification 


Unique sRNAs (%) 


Total sRNAs (%) 


Total sRNAs 


2917391 (100%) 


24449036 (100%) 


Common reads 3 


126807 (4.35%) 


17718517 (72.47%) 


A lumbricoides specific 6 


1711536 (58.67%) 


4379273 (17.91%) 


A suum specific c 


1079048 (36.99%) 


2351246 (9.62%) 



Note: a reads that were shared by the two taxa; b reads found in A lumbricoides, 
but not found in A suum; c reads found in A suum, but not found in 
A lumbricoides. 
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Table 2 Comparison of miRNA profiles in Ascaris 
lumbricoides (Alu) and A. suum (Asu) 





Shared 


Asu-specific 


Alu-specific 


Asu-total 


Alu-total 


Novel 


11 


396 


92 


407 


103 


Known 


63 


24 


5 


87 


68 


Total 


74 


420 


97 


494 


171 



Results 

Profile differences in short RNAs between the two 
Ascaris taxa 

Deep sequencing yielded 18.29 and 11.72 million raw 
reads in adult of A. lumbricoides and A. suum, respect- 
ively, with 14.69 and 9.76 million high quality reads that 
longer than 18 nt. Length distribution analysis showed 
that the reads of A. lumbricoides and A. suum were fo- 
cused on 21-23 nt. Among the clean reads, 46.11% and 
39.97% were identified as non-coding RNA (ncRNA) in A. 
lumbricoides and A. suum respectively, including tRNA, 
rRNA, snRNA and snoRNA, which were at near the same 
level. Among the high quality reads, 72.47% of the total 
reads were shared by the two parasites (Table 1), while 
1.71 and 1.08 million was A. lumbricoides and A. suum 
specific. 

miRNA profiles of the two taxa 

By mapping onto the A. suum genome, we obtained 171 
and 494 miRNA candidates, with precursors meeting the 
criteria listed in the method and having standard stem- 



loop structures (Table 2, Additional file 1: Table SI). 
Among the miRNA candidates, 68 and 87 miRNAs were 
matched with the known A. suum miRNAs deposited in 
the miRBase database, and 63 of them were shared. The 
left miRNA candidates were marked as novel miRNAs, 
and among which only 11 miRNAs were shared, there- 
fore, 92 and 396 miRNAs were A. lumbricoides and A. 
suum specific. Totally, there were 74 miRNAs shared by 
the two nematodes, including 63 known miRNA and 11 
novel ones; and there were 97 and 420 miRNAs being A. 
lumbricoides and A. suum specific. 

Target prediction and the functional prediction of the 
predicted targets 

A total of 57,359 mRNA and EST items of Ascaris de- 
posited in NCBI were downloaded and used for target 
prediction. Under the stringent matching criteria, it was 
found that the target numbers of both nematodes 
ranged from one to thousands. For A. lumbricoides, the 
target number ranged from one (Alu-miR-novel-012-3p, 
Alu-miR-novel-066-5p, Alu-miR-novel-102-3p) to 3,562 
(Alu-miR-novel-063-3p), 255 in average. For A. suum, 
the target number ranged between one (asu-miR-novel- 
027-3p, asu-miR-novel-044-5p, and asu-miR-novel-068- 
5p) and 3,343 (asu-miR-novel-039-3p), 168 in average. 

Functional prediction revealed a significant set of targets 
which are related to ovarian message protein in A. lumbri- 
coides (n = 9) and A. suum (n = 10), a set of vitellogenin (5 
in A. lumbricoides and 18 in A. suum), and a set of 



Alu_specif.txt Asu_specif.txt 




Cellular Component Molecular Function Biological Process 



Figure 1 Enrichment analysis of the functions of strain-specific miRNAs targets of Ascaris suum and A. lumbricoides. The horizontal axis: 
Gene Ontology analysis to the targets, including cellular component, molecular function, and biological process. The vertical axis: percentage of 
genes in total targets. 
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Figure 2 The main blast result and standard stem-loop 
structure of precursors of one of the representative miRNAs. 

Upper: plain text sequence of precursor of asu-miR-share-031 and 

the main blast result with mature miRNA matched; Down: stem-loop 

structure of the precursor. Detailed blast information was shown in 

(Additional file 2: Figure S1). 
\ / 



chondroitin proteoglycan (5 in A lumbricoides and 7 in A. 
suum). The same phenomenon is found in NADH de- 
hydrogenase (n = 7 in A. suum, n = 3 in A. lumbricoides), 
40s ribosomal protein (n = 14 in A. suum, n = 2 in A. lum- 
bricoides), and 60s ribosomal protein (n= 11 in A. suum, 
n = 7 in A. lumbricoides). As a distinguished character, 
there were more targets related to movement in A. suum, 
such as actin (n = 5) and tubulin (n = 6), while such targets 
were very few in A. lumbricoides. 

Gene Ontology (GO) analysis of the targets of miR- 
NAs resulted in 3 respect outcomes, including cellular 
component, molecular function and biological process. 
Enrichment analysis showed that percentages of most 
targets functions of A. suum and A. lumbricoides were 
similar, except a few items (Figure 1), which indicated a 



very closely metabolism pattern of the two parasites. 
However, for cell component, miRNA targets of A. suum 
had an extra cellular component part named as "synapse 
part"; for molecular function, one function named "elec- 
tron carrier" was only found in targets of A. suum. In 
addition, three biological processes named biological ad- 
hesion, death, and immune system process were only 
found in the miRNA targets of A. suum. 

Quantification of A. lumbricoides and A. suum miRNAs 

Three A. suum and A. lumbricoides specific novel miR- 
NAs, including asu-miR-novel-383, asu-miR-novel-097 
and asu-miR-novel-031 in A. suum, and alu-miR-novel- 
053, alu-miR-novel-021, and alu-miR-novel-064 in A. 
lumbricoides, were representatively selected for quanti- 
fication using modified stem-loop qRT-PCR (Figure 2). 
These 6 miRNAs had the lowest A G energy, and/or 
more matched variants, and had mature miRNAs at 
both of the arms of the stem-loop precursors. The rep- 
resentative blast and standard stem-loop structure as 
asu-miR-novel-031 was shown in Figure 3. Detailed 
blast information with variants of others was shown in 
Additional file 2. 

Both of the mature miRNAs of the 6 selected miRNAs 
can be successfully detected, except the alu-miR-novel- 
64-5p and asu-miR-novel-383-3p. Two miRNAs named 
asu-miR-novel-031 and alu-miR-novel-53 had significant 
expression differences at 5p and 3p. For alu-miR-novel- 
53, the expression levels were 1.00 ± 0.12 at 5p, and it 
was 1.68 ± 0.17 at 3p. And for asu-miR-novel-031, the 
relative expression level was 1.85 ± 0.24 at 5p, while it 
was 1.56 ± 0.30 at 3p. The alu-miR-novel-64-5p and asu- 
miR-novel-383-3p were not successfully amplified, which 
could be attributed to the very low expression levels, 
despite that the primers and reaction mixture were 
modified several times. 



Asu-miR-novel-031 gi | 32025639 1 1 gb |AEUI0 106 1677. 1 1 : 710 : 785:+ 76 (nt) -60. 00 (kcal/mol) 
GUUCAAGGAAAGGAUUGUGCAGACGUGAAAGCGUUCAACUAACGCUUUCACGUCUGCACAAUCCUUUCCUUGAUCU A sum iR-nove 1-0 3 1 348 
..((((((((((((((((((((((((((((((((( )))))))))))))))))))))))))))))))))... 

**********AGGAUUGUGCAGACGUGAMG Asu-miR-novel-031-5p 322 

********************************* Asu-miR-novel-031-3p 23 



Asu-miR-novel-03 1 -5 p 




Figure 3 Quantification of representative miRNAs of Ascaris lumbricoides and A. suum. Three novel miRNAs of each nematode were 
detected. Two miRNAs named asu-miR-novel-031 and alu-miR-novel-53 had significant expression differences of 5p and 3p mature miRNAs. 
The alu-miR-novel-64-5p and asu-miR-novel-383-3p were not successfully amplified, which might due to their low expression levels. 
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Discussion 

A previous study revealed that A. lumbricoides and A. 
suum have identical 5.8S and ITS -2 rDNA sequences 
[4]. There were only six (1.3%) nucleotide differences in 
ITS-L All these genetic evidences supported the pro- 
posal that A. lumbricoides and A. suum represent a sin- 
gle species [4,5]. 

Due to their specific expression in different organisms, 
tissues and cells, miRNAs may provide potential novel 
biomarkers [27-29]. Our previous study indicated that 
different T. gondii genotypes have their unique miRNA 
profiles [30]. In the present study, the two nematode 
taxa shared very high percentage of total reads (72.47%), 
while having very small percentage of unique reads 
(4.35%), and the total percentage of ncRNA was at simi- 
lar level (46.11% and 39.97% in A. lumbricoides and A. 
suum, respectively), which indicated huge redundance of 
some ncRNA of the two nematodes. When parasitic en- 
vironments change, the gene expression of parasites will 
be correspondingly changed, and regulators such as 
miRNAs will also be changed, therefore the death and 
new creation of miRNAs are very common [36]. For the 
different parasitic environment of pigs and humans, adap- 
tive modifications should have occurred. So although the 
sRNA/shared-miRNAs for fundamental metabolism were 
similar in the two taxa, there still are some specific miR- 
NAs in each nematode taxon. We obtained 171 and 494 
miRNA candidates as expression profiles, with 74 miR- 
NAs being shared by the two nematode taxa, and 97 and 
420 miRNAs being A. lumbricoides and A. suum specific. 
The same precursors of a miRNA can be found in the 
genome of A. lumbricoides or A. suum, however, we ob- 
tained fewer miRNAs from A. lumbricoides than from A. 
suum, which indicated that some miRNAs in A. lumbri- 
coides were not expressed. This difference might be re- 
sulted from the different parasitic living environment in 
pigs and humans. 

Function prediction and enrichment analysis showed 
that targets of the miRNAs of the two parasites had 
similar metabolism patterns, including cellular compo- 
nent, molecular function and biological process, with 
some specific differences, such as different target num- 
ber, specific functions (NADH dehydrogenase and elec- 
tron carrier functions), etc. 

Conclusion 

The present study characterized comparatively the miR- 
NAs of adult A. lumbricoides and A. suum, and the 
findings support the recent proposal that A. lumbri- 
coides and A. suum represent a single species [5]. Due 
to the fast evolution nature of miRNAs and the differ- 
ent parasitic living conditions of humans and pigs, the 
phenomenon above might indicate a fast evolution of 
miRNAs of Ascaris. 
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Additional file 1: Table SI. Shared and strain-specific miRNA with their 
predicated target of Ascaris lumbricoides and A. suum. 
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